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1 Introduction 



We shall deal with Feynman integrals corresponding to the so-called "sunset" diagram. 
This is a two-loop, self-energy-type diagram involving three propagators. In what follows, 
we shall need different integer powers z/j of the corresponding denominators. The scalar 
integrals corresponding to this diagram with an external momentum k are denoted as 



L(n;z/i,z/2,z/3) = J J 



d> (Tq 



[{p - g)2 - mlY' [g2 - [{k - pf - mlY' 



where n = 4 — 2e is the space-time dimension in the framework of dimensional regulariza- 
tion The sunset diagram possesses a three-particle threshold at fc^ = (mi +m2 + m^Y 
and (in general) three pseudothresholds, at k"^ = (mi -|- m2 — ms)^, (mi — m2 + m^)"^ and 
(-mi + m2 + ma)^. 

There are several reasons why analytical evaluation of such diagrams and, in particu- 
lar, its threshold and peudothreshold values, is important: 

(i) They are needed for calculation of some realistic radiative corrections in the Standard 
Model and its extensions (in particular, in the Higgs sector). The sunset-type integral is 
a part of the basis for any two-loop two-point calculation (see e.g. in 0). 

(ii) Although this is one of the basic two-loop-order diagrams in the Quantum Field The- 
ory, analytic results for general values of the external momentum k and the masses m^ 
are not available]^, at least in terms of special functions like (generalized) polylogarithms. 
Moreover, there are some arguments P] that the results for such diagrams cannot be ex- 
pressed in terms of polylogartihms, with a possible exception of special values of k"^. 

(iii) This is the simplest example of a diagram involving a three-particle cut with all the 
three particles being massive. Exact results for this diagram would be useful for under- 
standing the analytic structure of three-particle cuts. Note that the four- dimensional case 
is much more complicated than the three-dimensional one considered in 

(iv) The threshold values of such diagrams are needed for constructing analytic approxi- 
mations of the behaviour near the threshold, cf. e.g. refs. |^ ^ |TU|. 

(v) In some cases (for example, when two masses are equal, mi = m2), the pseudothreshold 
may coincide with the on-shell limit fc^ = m| which is relevant for the on-shell calculations 
(see e.g. in refs. 0). 

(vi) Some numerical approaches (in particular, the one described in [|l^) involve the 
sunset-like integrals as "kernels" of integral representations for more general two-loop 
graphs. It should be noted that some integral representations which can be used for nu- 
merical calculation of two- loop self-energy diagrams can also be found in refs. [13, IH]. 



(vii) The threshold value corresponds to an infinite sum related to the small momentum 
expansion []TB|, |TB|, or the large momentum expansion ITU], taken at its border of con- 



vergence. The closed form for the coefficients of such expansions (in terms of generalized 



hypergeometric functions) was given in ref. ||2^. Therefore, as a by-product of the calcu- 



lation one could get some summation formulae for complicated hypergeometric functions. 
Furthermore, since the sunset diagram may be represented in terms of a one- dimensional 
integral involving four Bessel-type functions |^ the analytic results are applicable 
to those integrals, too. 

^The results for some other two-loop self-energy diagrams can be found in refs. ^, ^, ||, ||]. 
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2 Approach to the calculation 

The standard Feynman parametric representation for the integral (|l]) reads 



1 1 1 



L[n] Z/i, Z/2, Z^sj = 1 TT 



r(z/i) r(z/2) r(z/3) 7 7 7 (aia2 + ttias + a2a3)^''/^ ^''^ 



[010203^^ — (tti«2 + ctias + a;2«3)(aimf + 0:2^2 + asm!)] 



(2) 



Using tricks similar to those described in we can get rid of the first denominator in 
the integrand of eq. . Namely, let us use exactly the same transformation of a- variables 
(first inverting and then rescaling) as in eq. (13) of ||22|| , 

= («D"\ = ^(a"i, a"2, «"3), ^(«"i, a"2, ^"3) = { ^'^''^ ■ (3) 

a I + a 2 + ot 3 

Suppressing the primes, we arrive at the following modified representation: 

T I \ •2-2n n ^(i^l \ V2 ^ V'i - u) 

L(n;z/i,z/2,z/3) = 1^ ^" vr" ' 



r(z/i) r(z/2) r(z/3) 

111 i^2+!^3-n/2-l j/i+i'3-n/2-l i/i+i/2-n/2-l 



^ . . . Qi ' ^2 ndaj (5(Ea» - 1) 

Ill r r n O O 2l^'^i~"- ' 







The equivalence of the representations @ and (^) can also be established by comparing 
the corresponding triple Mellin-Barnes contour integrals. 

Furthermore, using Cheng- Wu theorem (see also in , Appendix B) and rescal- 



ing the variables, one can transform the representation (^ into 

.2-2S1., r(l/i + U2 + iy3-n) fyj ^n/2-u. 



T{ui) T{u2) T{u3) 



U2+U3-n/2-l ^ui+us-n/2-1 



(mi^+m2?7+m3)3"/2"^'^' [{rai^+ra2r]+ra3){niir]+m2i+ra3^7]) — k'^^rjf^'^''^ 





The cubic form in the denominator has the following representations at the pseu- 
dothresholds and at the threshold: 

{miS, + 11127] + m3)(mi?7 + m2^ + nisirf) - k'^^i] 

= 77127713^(1 - 77)^ + 7711771377(1 - + 771i77l2(^ - 7/)^ + (^(771i + 7712 + "^3)^ " 

77l27H3^(l — 7/)^ + 7111711377(1 + + 711i71l2(^ + 7/)^, /e^ = (— THi + 7112 + ^3)^ 



71127113^(1 + 77)^ + 711i71l37/(l — + 711i71l2(^ + 7/)^, /e^ = (tIIi — 7II2 + Tlls)^ 
71127113^(1 + 77)^ + 711i71l37/(l + + 711i71l2(^ — 7/)^, /e^ = (tIIi + 7II2 — Tlls)^ 

71127113^(1 — 77)^ + 7111711377(1 — + mim2{^ — riY-, k"^ = {nii + 7112 + 7113)^ 



(6) 



In particular, it can be seen that this cubic form is positive semidefinite at the threshold 
and pseudothresholds. 
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Using eq. @, it is easy to derive the following decomposition: 

L(4 - 2e- 1, 1, 1) = {-k\-^L{Q - 2e; 2, 2, 2) 

+mlL{A - 2e; 2, 1, 1) + m^L(4 - 2e; 1, 2, 1) + mlL{4: - 2e; 1, 1, 2) | . (7) 



This decomposition is similar to one used in |]T3[. Taking into account that 



L(6 - 26; 2, 2, 2) = -vr'^Iv(4 - 2^; 1, 1, 1), (8) 

eq. d^) just reflects the fact that the mass-squared dimension of the integral L(4— 2^; 1, 1, 1) 
is {l-2e). 

The calculation of each of the integrals on the r.h.s. of is simpler than direct 
calculation of L(4 — 2e; 1, 1, 1). Using the representation (^), we see that the integrals 
L(6 — 2e; 2, 2, 2) and L(4 — 2e; 1, 1, 2) are proportional to 

+ ?Ti2?7 + rris) [cubic tormj 







where r = 3 for L(6 - 2e; 2, 2, 2) and r = 2 for L(4 - 2e; 1, 1, 2). 

When r = 3, the double integral in is convergent, and what we need is just to 
expand the integrand in e, keeping the terms of order e (since we have got a singular 
factor r(2e) in front of the integral). In this way, we get 

f f dq dri 

r(2e) / / {1 + eln^ + eln?7 + 3eln(mi^ + m2?7 + ms) 

-2eln[cubic form]} + (10) 







When r = 2, the integral (^ develops a (l/e) singularity as ^,?7 oo. One can 
subtract from (^ a simpler integral with the same asymptotic behaviour as ^,?7 ^ oo. 



d^dr] ^^7]^ 



2e ' 



which can be calculated in terms of F functions. The difference of (^) and (11) is conver- 
gent and yields 

d^dr^ ^ /(mie + m2r^)m3er/\ ^^2) 







{vTLx^^ + m2ri -|- ma) 2 y [cubic form] y 



When we consider representation (|) for L(4 — 25; 2, 1, 1) and L(4 — 2e\ 1, 2, 1), it does 
not diverge as ^,?7 — oo. Instead, it does diverge as ^ — or 77 — * 0. However, these 
integrals can be reduced to L(4 — 2e\ 1, 1, 2) just by a permutation of mi, m2, ms and thus 
calculated via subtractions similar to (|Tl|). 



4 



The most complicated point in calculating the parametric integrals (0) and ([1^) is 
how to deal with the contributions involving In [cubic form]. One of possible ways is to 
substitute the variables, 

^ = 7 (m3 + m2p), r] = {ms-mip), (13) 

mi + m2 mi + m2 

and then integrate over p between (—7713/7712) and (m^/mi). Depending on a region 
(and whether a pseudothreshold or the threshold is considered), this integral yields an 
arctangent or a hyperbolic arctangent (the latter can be presented as a logarithm) of an 
argument involving square roots. The most labour-consuming part was to calculate the 
remaining integral over a, which was done (for the cases considered below) by using some 
tricky trigonometric substitutions. 



3 General results at the pseudothreshold 

We shall consider the pseudothreshold at A;^ = (r77i + 7^2 — 7/23)^. The results for the other 
two pseudothresholds, at k"^ = {mi — m2 + 7/73)^ and A;^ = (— r77i + m2 + m^Y can be 
obtained by permutation of the indices 1,2,3. 

It is convenient to introduce the following dilogarithmic functions^: 

T-{z) = U2i-z) -U2{-l/z) + \nz\n[{l + zf/z) 

= 2U2{-z) + Itt^ + 2\nz\n{l + z) - ^In^ z 

= 2U2{l/{l + z))-l7T^ + ln\l + z)-l\n^z, (14) 



T+(z) 



U2 (1 - 1/z) 
2U2 (z) 



3^ 



Li2 (1 - z) 
f 21nzln(l - z) 



Un'z 



2Li2(l/(l-z)) + iln2 



(15) 



which are antisymmetric under inversion, T^(l/z) = — T^(z). In particular, T^(l) = 
and T^{mj/mi) = —T^{mi/mj). Note that the functions T^[z) continue each other to 



the region of negative z's, namely T 



±1 



T^{-z). 



Let us list the results for the integrals involved on the r.h.s. of eq. (0): 



L(6-2e;2,2,2) 



TT 



l+e) 



1 m 

'4e 



m 



mi + m2 - 7773 



9 

^ + 



m\ + 7172 + 777| 

2(7/21 + — nT'sY 



7r7i7r727r73 



{mi+m2 — m^)^ 
222 

(7771+7772-7773)4 



7^71 — 7772 



7773 



m^ 



-mo 



m^mo 



In 


mi 




7772 




_ /777i 


T 






\7772 



7772 + 7773 _ 7772 

m — 



777i+7r73 , 777i' 

m — 



7r7i 



7773 



7772 



7773 



+ 



7r7n 



■7r7 



2 2 

m2m^ 



^ T+{ — 
7773 



7r7o 



■7r7 



m^m 



^Up to In^ z terms, these functions are similar to those used in refs. iQ, |ll|, [l^ . 



1 r+ {"^ 

mi 



5 



rrin 



2ml 



2ml 



+ 



m\ 2m^ 2ml 



^ m2 3 ^ 



ml 2ml 2mf 



In 



2 mg 

mi 



■ + 0(e),(16) 



L(4-2£;l,l,2) 
1 + ^ 



TT 



4- 2s 



r2(i + £) 



-4e 



1 

2? 



Vma/ 2 ms 2 ms 



Vms 



mi+m2 — m3 
1 

(mi+m2— ms)^ 



mi ma 
mi m h m2 In — 

ms ms 



,m 



^, , mi 7r^\ 2m^ — m? n mo 2m?— m^ , r, mi 
In^— + ^ + ' M n^^+ ^ ^ In^^ 
, m2 o / 2 ms 2 ms 



+ 0(e),(17) 



L(4-2£;2,l,l) 
f 1 



TT 



A;^=(mi+m2— ma)^ 
2 



o + 1 T- 

2^2 6 



mi\ 
m2/ 



1-e 
/miN 1 
Vms/ 2 



m 



-4£ 



In^ — + - In^ — 



m2 



ms 



+ 



mi+m2 — ms 
1 

{mi+m2-m-i)' 



mi ms 
m2 In h ms In — 

m2 mi 



m. 



2ml — ml A 2 "^1 



In"— + — - 
m2 3 y 



ml+ml 2 "^2 2m|— m2 2 "^i 



In" — + 

ms 



In" — 

ms 



■ + 0(e).(18) 



The result for L(4 — 2e; 1,2, 1) can be obtained from 
Substituting eqs. ([T6|)-(p!8D into eq. we get 



by permutation mi m2. 



L(4-2£; 1,1,1) 
1 f 2 



4-2e 



r2(i+£) 



A;2 = (mi +m,2 — m,3 ) 



2^2 



2-4e+^2-4e 



(1_^)(1_25) 

+ ml~^') + —(mi + m2 - ms) (m\-^' + m\-^' - m\-^') 



7 1 TT^ 

+ -(mi + m2 - ms)^ + -{m\ + m\ + ml) — —{jn\ + m2 
o 2 



~{m\ — m2)T 



- / ^1 



(m2 - m3)T+ 
m2y \m^J 
2 I ^2 ^ ^2 



m2\ _ 



mo 



m^)T^ 



ms 
mi 



mj + m| j^2 ^ _^ "^2 + "^3 j^2 m2 _^ + mf ,^2 "^3 



+ 



mim2ms 
mi+m2— ms 



m2 
mi — m2 

ms 



In 



2 

mi 

m2 



ms 



mi 



2 

ms 



In"^ 

mi 



m2 + ms m2 mi + ms mi 
In In — 



m2 



ms 



6 



2 2 2 



m\ — m\ ^_ /nil 



ml- ml 

' 2 2 

m2m^ 



m2 




m\-m\ ^ 

2 2 



2ml 



1m\ 



In 



2 ^ 



"^3 

mi. 



-0(5).(19) 



4 General results at the threshold 



In addition to the notation (14), we shall also introduce the angles ^i, ^2 and 6'3 via 



Qn = arctan m^ 



Imi + m2 + ms 
mim2m^ 



(20) 



Note that 9i + 62 + 9^ = vr. Therefore, they can be undestood as the angles of a certain 
triangle (cf. e.g. in [p5| ). 

The results for the threshold values of the integrals occurring on the r.h.s. of eq. (0) 

are 



L(6-2e;2,2,2^ 



TT 



X ■ 



1 m 

'4e 



1-4^ + ml-^' + mj-^' 
mi + m2 + m^ 



9 



+ 



m^ 



m% 



8 2(mi + m2 + m^) 



mim2m3 



mi —m2 



{mi+m2+m3y 

222 

(mi + m2+m3)4 



+ 



mi — m3 
m 1 



m2 m^ — mi 
m 1 



In^' 



2 2 
m^ml 



-Ait 




In^^ 



-^ + -% + -% I +47r ( — + — + — 

mf 777.2 '^S / ^'^l ^2 



mim2m^ 



7773 
777i 



-0(e),(21) 



L(4-2£;l,l,2: 



-2i^ + ^+ 3 



fc^=(mi+m2+m3)^ 

47r2 



TT 



777i 



4-2£ 



r^(i + e) 

l-e 



777. 



-4£ 



7773 



7772 



7773 



7771 +7772 + 7773 
1 

(7771+7772 + 7773)2 



7771 7772 

777i m h 7772 lH 

77I3 777,3, 



1 2 "^1 1 2 ^2 

- m 1 — m — 

2 7773 2 7773 



777^T- f ^ 



7773 



7 



mf + 1712 Y 2 "^1 _|_ — ^2 J 2 ^1 _|_ — mf ^ 2 "^2 



1712 



7713 



77I3 



47r2 



(^772^+7712^ —Att (^mg^i +777^^2) —47r-y/mim2m3 (?77i +7772 +7773) 



■ + 0(5).(22) 



The results for L(4 — 26; 2, 1, 1) and L(4 — 2e; 1, 2, 1) can be obtained from p^ ) by per- 
mutations 7773 ^ 777i and 7773 777.2, respectively. Using eq. we get 



L(4-2e; 1,1,1 
x{-^ H 



,4-2e 



THl + e) 



TT 



'l-e){l-2e) 



2-4.+^2-4.^^2-4e 



47r2 



+ ^("^1 + "^2 + "^3) ("^i' 



+ 777^-^^ + 777^-^^ 



+ -(7771+7772+7773) +I2+ 3 



(777^ + 7772 + "^3) - 477 (777^6'i + 77726^2 + 77736*: 



777 



777n 



ml)T' 



/7772 
V7773 



777o 



777?)T- 



7773 
777i 



mi+mi 2 ^1 7770+7773 2 "^2 mi+mi n "^3 / ^ \ 

In 1 In 1 In 477-1/771177727/73(7771+7712+7773) 

2 7772 2 7773 2 777i * 



+ 



77117/727/73 

777i + 7^72 + 7773 

2 2 2 
7/7^7/727773 

(7771+7772+7773)^ 



7/71 — 7772 , 777 i 7772 — 7/73 7/72 ^-3 — T^l , TTls' 

m 1 In 1 In — 



7773 

2 2 
777 J— 7772 

2 2 
777 i7/72 



7772 ^1 
/777i 

\7772 

1 



7773 



7772 



T 



+ 



m\-77ll ^_ ^7772 
77727773 



V7773 



+ 



777i 

ml-77il 

2 2 
7/7i 7773 



"^3 
777i 



2777i 27^72 



In^^ 

7772 



771 



27772 277i 



In^^ 

7773 



+ 



^7772 27^73 2777i^ 
7/7f 7/72 '^i 



. 7771 47r V 1 1 1 

3 Y7/7i 7772 ^3, 



-47r I — + — + — f I +47r ( — + — + — 



7773 

1 1 1 
— + — + — 

777 i 7772 ^3 



/ 7771 + 7772 + 7773 
7/7i 7/72 7773 



'+0{e){23) 



The presented expressions explicitly obey all required symmetries. In particular, the 
result ( p2|) is symmetric under 7771 ^ 7772, whereas eqs. (^l]) and (|2^) are totally symmetric 

in 777i, 7/72, ""T-S- 



5 Some special cases 

When two of the masses are equal (say, 7771 = 7/72 = 7/7, 7/73 = M) the presented general 
expressions can be simplified. Let us consider, as an example, the integral L(4 — 25; 1, 1, 1) 
and introduce a dimensionless variable /i = M/7/7. 

At two different pseudothresholds, fc^ = {2m— = (2— /i)^777^ and k'^ = = /i^7r7^, 
we get the following results: 



L(4- 25; 1,1,1) 



= r (1 + g) ^2-4. 

mi=m2=m, m,z=M=^in (1 — ^)(1 — 2^) 

k'^ = {2m-Mf = {2-n)'^mP' 
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x|-^ + L_J^ + ^ln,+ -(2 + ,^) + -(2-,f + ^^,ln, 



(2-/i)^ 



Li2(l-/i) + 



12 



(24) 



L(4-2£; 1,1,1) 



77 



4-2£ 



r^d + e) 



mi=m2=m, mg=M=iMm 



:i-e)(l-2e) 
11 



2-4e 



^"^^ +^ + -/i^ln/i + l + - (2 + /i^) In /i - 2/i^ In^ /i 



2^2 4e e' 



:i-/i 



2\2 



Li2(l-/x2)-y }+0(^). (25) 



Note that the latter expression is relevant for the on-shell calculations and corresponds 
to the results presented in refs. [|ll|, [1^ . 

At the threshold, when k"^ = (2m + Mf = {2 + n) ^m^, we get 



L(4-2£; 1,1,1) 

f 2 + (2 + ^)2 2/i2 



mi=m2=m, m3=M=iim 
fc2 = (2m+M)2 = (2+At)2m2 



71 



4~2e 



;i-5)(i-25) 



m 



2-4e 



X 



+ 



2^2 ■ 4£ 

,2 



-ln/i+^(2 + /i2) + I(2 + /i)2 



6 + 2/i + yU^ , o 21 2 2 

-/iln/i — 2yU In fJ' — r:77^(3 + 2/i + yU jTT 



2 + /i 

4(l+^)2(l-/x)(3+/i) 
(2 + /^)^ 



Li2 



(2 + /i)3/2 
1 \ 1, 2. ^ /I 

+ -ln (l+/i)+7rarccos 



l+/i/ 2 



l+/i/ 12 



■ + 0(5).(26) 



For /U = 1, i.e. when all masses are equal, the results (|^) and (^) give the same as 
eq. (30) of the paper |^ (see also in |^), 



4-2. r\l + e) ^2-4er_ 3 ^ 1 ^ 19- 



L(4 - 2e: 1, 1, 1)1 = vr^ — ■ — -m , , , 

^ ' ' ' ^' mi=m2=m3=m (l-£)(l-2£) I 2^2 Ae 



(27) 



whereas eq. (p6D yields the threshold value of the sunset integral with equal masses 
L(4-2£;l,l,l)| 



TT 



4-2. r2(l + e 



m\=m2='mz=m 



2 3 9 75 Svr 1,, , 



-m 



We have also reproduced the same result (E3) by using the dispersion technique. 



(2^ 



6 Conclusions 

We have obtained analytic results for the threshold and pseudothreshold values of the 
sunset diagram with arbitrary masses. The results are expressed in terms of dilogarithms 
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of ratios of the masses, all other functions being elementary. The results are explicitly 
symmetric with respect to all symmetries they have to obey. Certain checks, including 
numerical ones, confirm that the results are self-consistent. Note that most of the terms 
in eqs. (^)- (|23|) look just as naive analytic continuation of the results (p!^)-([T9|), if one 
changes into (—7713). The non-trivial terms are those involving vr, including the 6i 
terms. 

We note that the complete (unsubtracted) results, involving the and l/e ultra- 
violet poles, are simpler than the finite function obtained by subtraction of two terms of 
the expansion in fc^ near k"^ = (cf. e.g. in |]T3[). In the latter case, we would also get 
dilogarithms (or Clausen functions) corresponding to the vacuum diagrams |2^, |16|. In 
particular, for the equal- mass case ( ^8]) the Cl2(7r/3) contribution is missing in eq. (|28|) . 

As an extension of the presented approach, two-loop self-energy diagrams with four 
and five propagators may also be considered. For example, the exact result for the pseu- 
dothreshold value of the "master" diagram with equal masses was presented in p^. Al- 



though the results obtained may be used for all purposes listed in the introduction, we 
would like to emphasize that their application to constructing expansions near the thresh- 
olds is one of the most important points. 
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